An optimal topological spin pump 
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We study the recently introduced Z2 pump consisting of a family of one-dimensional bulk insu- 
lators with time reversal restriction on the pumping cycle. We find that the scattering matrices of 
these pumps are dichotomized by a topological index. We show that the class of pumps character- 
ized by a nontrivial topological index allows, in contrast to its topologically trivial counterpart, for 
the noiseless pumping of quantized spin, even in the absence of spin conservation. This distinction 
sheds light on the Z2 classification of two-dimensional time reversal invariant insulators. 
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Introduction: — The idea to pump spin through a meso- 
scopic device at zero bias by cyclic variation of two sys- 
tem parameters is very appealing due to its promise of 
precise and reversible flow control [1 . Optimal pumps 
which operate noiseless, transferring quantized spin in 
a cycle are particularly relevant for potential applica- 
tions p5. In this paper, we propose the concept of an 
optimal topological spin pump, composed of a bulk in- 
sulator. We show that the ability to pump quantized 
noiseless spin is a hallmark of the nontrivial topological 
invariant characterizing a quantum spin Hall system. 

The discovery of the quantum Hall effect introduced an 
interesting new way to classify different states of matter. 
Unlike more familiar phases, a quantum Hall state does 
not break any symmetry and cannot be described by a lo- 
cal order parameter. Rather, it differs from a regular two- 
dimensional insulator by topological invariants, known as 
Chern numbers, which reflect the global structure of its 
ground state wave function [31 [4]. This topological clas- 
sification of states of matter found recently an exciting 
extension to two-dimensional time-reversal invariant bulk 
insulators [5]-[7]. In accordance with the Chern numbers 
used to classify quantum Hall systems [6 , these systems 
can be characterized by a Z2 topological index, based on 
the properties of their bulk ground state. The group of 
insulators described by a nontrivial Z2 index are known 
as quantum spin Hall systems. 

In analogy to the quantum Hall state, it is possible to 
gain insight into this topological classification by study- 
ing a pump formed by placing the two-dimensional sys- 
tem on a cylinder threaded by a magnetic flux |9]. 
The resulting system can be mapped onto a set of one- 
dimensional time dependent Hamiltonians by identifying 
the magnetic flux with time. This mapping defines a Z2 
pump [7 . Nonetheless, in the absence of spin rotation in- 
variance, the physical meaning of the pumped Z2 charge 
remains elusive. In this letter we study the scattering 
matrix of the open one-dimensional pump constructed 
in this way. We find that mapping the two-dimensional 
insulating system onto a one-dimensional pump estab- 
lishes a Z2 classification of the scattering matrix. The 



resulting one-dimensional pumps are dichotomized. The 
family of scattering matrices belonging to the topologi- 
cal nontrivial class allows, in contrast to its topological 
trivial counterpart, for the noiseless pumping of a quan- 
tized spin, even in the absence of spin conservation. We 
illustrate these ideas by two examples. 

'L2-index: — We consider a family of one-dimensional 
Hamiltonians of non-interacting electrons with a bulk en- 
ergy gap, that depend continuously on a cyclic pumping 
parameter t, H{t + T) = H{t)^ and satisfy: 



H{-t) 



GyH- 



\t)(Jy, 



(1) 



where a^ are the Pauli matrices. These systems can be 
viewed as a mapping of the set of two-dimensional time- 
reversal invariant insulators placed on a cylinder, where 
t corresponds to a magnetic flux threading the cylin- 
der. Indeed, Eq. (IT]) implies time-reversal invariance of 
the corresponding two-dimensional system, as is evident 
upon identification ikx^t) -^ (kx^ky). In the context 
of pumping, Eq. (IT]) ensures the existence of two time- 
reversal invariant moments (TRIM) ti = and ^2 = T/2, 
at which H{ti) = (TyH^ {ti)(jy where i = 1,2. Upon cou- 
pling the one-dimensional system to two single channel 
leads, the open system is described in terms of a time- 
dependent 4x4 unitary scattering matrix. Provided the 
system exceeds the attenuation length associated with 
the bulk energy gap, the transmission vanishes, and the 
scattering matrix is block diagonal. Each block, Tq,, is 
a unitary 2x2 reflection matrix in spin space, where 
the index a = L^R refers to the left and the right lead, 
respectively. The average spin injected into lead a dur- 
ing a cycle can be expressed in terms of the spin cur- 
rent, Sa = Imti (^[dra/dt]arl^) [lOj, see ^ below. It 
can be readily verified that Sa is invariant under gauge- 
transformations Tq, -^ e^'^^^^^fo;. This allows us to re- 
strict our analysis to the U{2)/U{1) :^ SU{2) particle- 
hole symmetric part of r^ which we denote by r^ • 

The pumping cycle defined in Eq. ([T]) provides a map- 
ping from the periodic time, corresponding to a one- 
dimensional loop, S^, to the SU{2) space of particle- 
hole symmetric reflection matrices, r^- Since the SU{2) 



group is topological equivalent to the three-sphere, 5^, 
one may conclude that a general pumping cycle does not 
distinguish between different topological classes, since all 
closed contours on S^ can be contracted to a single point. 
(The first homotopy class of S^ is zero, ni(S'^) = 0). 
However, Eq.(fT]) restricts the mapping by associating 
S{T/2 + t) = GyS^{T/2 - t)Gy. In particular, at the 
two TRIM, the SU(2) part of the reflection matrix is 
given by ra{ti) = il- This restricted mapping has two 
topological distinct classes, which may be categorized by 



f,(0)r„(r/2) = (-l)-^!. 



(2) 



Any loop f^ (t) on the three sphere characterized by i/ = 
can be contracted onto a single point, while paths with 
u = 1 are fixed by two distinct points at TRIM and can- 
not be contracted. Following (IT]), any additional points 
at which Tq, = — 1 away from the TRIM occur in pairs. 
We may, therefore, equivalently define v as the parity of 
Tq, = — 1 moments traversed in a cycle. Away from the 
TRIM, however, these points are not protected by sym- 
metry and can be removed by a perturbation larger than 
the level broadening introduced by coupling to leads, see 
Eq. ([5|. Hence, in the weak coupling limit, we can disre- 
gard any such accidental points. 

Topological spin pump: — The Z2 classification of the 
scattering matrix has a direct effect on the spin pumped 
during a cycle. We show below that the family of one- 
dimensional pumps belonging to the topological nontriv- 
ial class can, in contrast to its topologically trivial coun- 
terpart, operate as optimal pumps. To this end consider 
the effective Hamiltonian, H = Hl^Hr^ for the left and 
right edge states of the bulk insulator 



Ha{t) = fla{t) ^ ha{t) 



L,R. 



(3) 



Here we have employed the fact that left and right edge 
states decouple and the Hamiltonian is a sum of two 
independent 2x2 hermitian matrices in spin space, 
which can be parameterized by a unit matrix and the 
3 Pauli matrices, a. For transparency of the arguments 
we first consider //q, = 0, where the spectrum of (|3| 
is particle-hole symmetric, with eigenvalues ±\ha{t)\ = 
±/icK that vary during the course of a pumping cycle. 
Using the general expression for the scattering matrix 
5' = 1 + 2i7rWt (^H - iirWW'^y^ W, where W is a ma- 
trix that describes the coupling between the leads and 
the insulator [11 , and assuming that the leads couple 
equally well to up and down spins, so that Wa = i^^l, 
the reflection matrix of each block is given by 



m = m 



J(l){t)e(f){t)-a 



(4) 



where we have dropped the lead index a for brevity. Here 

Ht? - r^ . , , , 2r/i(t) 



cos (f>{t) 



/l(i)2+r2' 



sin(p{t) 



hity + r2 



(5) 



where F = it\w\^ is the level broadening due to coupling 
to the leads and we introduced the unit vector, 



e^{t) = Ht)/h{t), 



(6) 



which defines a time dependent rotation axis. We note 
that while the angle (j) depends on the level broadening, 
the vector e^ is exclusively determined by the effective 
Hamiltonian. Eq. ([5| shows that whenever the edge 
state crosses the Fermi level h{ti) = 0, the angle ^ = tt, 
and the reflection matrix is resonant f = — 1. 

Following Ref. [10] we express the spin injected into 
lead a during a cycle in terms of the reflection matrix 
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In general, the rotation axis e^ varies during the cycling 
process. As a result, one cannot identify a time inde- 
pendent axis along which spin pumped during a cycle is 
quantized. The situation, however, changes when the 
coupling to the leads is weak. We consider the limit 
1/T <C r <C E/^ where the level broadening F is small 
compared to the gap, yet large compared to the pump- 
ing rate, in order to allow for the pump to relax between 
cycles. In this limit the angle, (Isl), remains close to zero, 
(j){t) — 0(r/£^A), and changes rapidly to tt whenever 
a gapless edge state appears in the course of a cycle. 
The time duration of this transition can be estimated as 
^t - V I (dh(ti) I d£) - (r/£;A)T <C T. In the weak cou- 
pling limit, ^t can be made arbitrarily short in compar- 
ison to the time scale for variations of e^, so that e(/)(t) 
may be approximated by its value at the center of the 
resonance. The class of pumps with u = 1 cross a sin- 
gle resonance at the TRIM. Hence, in the weak coupling 
limit the reflection matrix of a topologically nontrivial 
pump describes a rotation around a fixed axis. 



r(t) 



i(f){t)e(j){ti)-a 
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(8) 



where ti is the TRIM at which the resonance occurs. As 
a result, the spin injected into lead a by the class of 
topologically nontrivial pumps is quantized 



S = he^{U){1^0{T/EA)) 



(9) 



where the quantization axis is determined by microscopic 
details of the system. Conversely, the class of pumps with 
u = either remains insulating during the entire cycle, 
or traverses two resonances at the TRIM. In the weak 
coupling limit, the former group can be approximated 
by a constant reflection matrix r{t) ^ 1, and thus does 
not pump spin. The second subgroup crosses two reso- 
nances during each pumping cycle. Each such resonance 
is associated with the value of the vector Cfj^iti). The two 
vectors 60(0) and e0(T/2), however, need not be aligned. 
As a result, one cannot identify a time independent spin 



direction that would lead to a quantized spin pumped 
trough an insulator with u = 0. 

These observations have a direct implication on the 
spin noise. Following Ref. [T2j [13], the variance of the 
spin pumped during a cycle, in a direction e^, is deter- 
mined by the time dependence of the vector ng - a = 
r(t)Ve* r{t) and vanishes for a constant n^ [13]. In 
the weak coupling limit, the reflection matrix of a topo- 
logically nontrivial pump, Eq. |8|, describes a rotation 
around a fixed axis, e0(ti). As a result, the vector n^^(t) 
with eq = e(f){ti) remains constant during the course of a 
pumping cycle. It follows that the topological nontrivial 
spin pump allows for the noiseless pumping of quantized 
spin. Conversely, a trivial insulator that crosses two reso- 
nances during a cycle cannot be parameterized by a time 
independent vector. Consequently, the trivial pump in- 
evitably operates with generation of finite noise. 

The above arguments generalize to finite chemical po- 
tential jiia- Once particle- hole symmetry is broken, the 
energy levels cross the chemical potential at two different 
moments, h{t±) = ±/i. From Eq. ^ it follows that these 
occur symmetrically around the TRIM, t± = ti± St. By 
diagonalizing the Hamiltonian at the crossing point, one 
can show that the width of each transition is determined 
by F, and in the weak coupling limit, is associated with a 
fixed vector e0(t±). We note that due to Eq. ([l]), e(/)(t+) 
and 60 (t_) are co-linear. Hence, r{t) describes a rota- 
tion around a fixed axis, and the spin pumped through a 
topologically nontrivial pump is noiseless and quantized 
in a direction which also depends on ja. 

Two examples: — To illustrate these ideas we next con- 
sider two examples which demonstrate the difference be- 
tween the two topological classes of pumps. To model a 
topologically nontrivial pump we consider the particle- 
hole symmetric Hamiltonian H = Hq -\- Vb -\- V^t -\- Vso 
[a [H] with 

n = &(t)E(-l)'<W.a^^./5 

i,a 
^so = Yl ^^so^a,/3 (4,aCi+l,/3 - 4+l,aCi,/3) • (10) 

Here sums run over TV sites and spin indices, and rst{t) = 
r°cos(27rt/T), b{t) = W sm{27rt/T). To simplify the il- 
lustration we choose ego = e^oX, such that r{t) depends 
only on two Pauli-matrices and can be represented as a 
point on the two-sphere. Finally, first and last sites are 
coupled to the leads via spin-independent hopping ele- 
ments Wa = '^al, while electrons inside the leads are 



described by a tight-binding model, Hq in (10). 

The transmission and refiection coefficients, Irl and 



f^, are found from the transfer matrix Trl by solving: 

^g^fcfio(A^+2)^^ ^ ^ 1 + Tl 



tRL [ jkRaiN+i)j — Trl 



-ikr 



eikLaf^ 



(11) 



where kR^L ^ ±7r/(2a) are the Fermi wave vectors in 
right and left lead at half filling, and the transfer ma- 
trix relates the states ('0Ar+2, V^at+i) = Trl (V^o, V^-i) , 



and is derived from Eq. (10). The resulting Irl is ex- 



ponentially suppressed over the length associated with 
the gap. To leading orders in Sr = Tst/ro, Sb = 6/ro, and 
Se = ego/To the refiection matrix is given by Eq. (pi with: 



0, 



/3i 



-Qfl 



.MTM W^^T^~M 



h 



V^^^M' 



(12) 



(ReAi + 25t) and pi 



Here we introduced ai 

(ImAi + 2(5e), with Ai = ^fW + 4(5r + i5ef. 

Fig. [l] shows the time evolution of the refiection ma- 
trix, r{t) = r(t)^ of model (10) with and without spin 



orbit coupling. The inset shows the time dependence of 
(j){t) and the y component of the vector e(/y{t)^ at finite 
spin orbit coupling and for different coupling strengths. 
At the TRIM, t = 0, T/2, the phase takes the value ^ = 
and (/) = TT, respectively. Consequently, r{t) belongs to 
the class of pumps with nontrivial Z2 index f = 1^ see Eq. 
([2|. The loops that represent the time evolution of r{t) 
are, therefore, fixed by the values at the two TRIM, and 
cannot be contracted to a single point by a continuous de- 
formation of the microscopic Hamiltonian that preserves 
condition ([l]). In the absence of spin orbit coupling, r{t) 
follows a geodesic, corresponding to the rotation around 
a fixed axis (black curve). At finite spin orbit coupling, 
the reflection matrices trace out curved loops, indicat- 
ing the absence of a fixed rotation axis. At moderately 
weak coupling, F/tq = 0.1 the refiection matrix follows a 
geodesic during most of the cycle (red curve), whose tan- 
gent is determined by the spin orbit vector ego- In the 
limit F/tq -^ the loops converge to a geodesic, implying 
that the pump works optimally. 

To model a topologically trivial pump we consider the 



Hamiltonian (10) with a double spatial periodicity of the 



time dependent parameters 

Vh = Kt) Yl (-l)''^«,/3(4n,aC2n,/3 + 4n+l,aC2n+l,/3) 



n,a,(3 



l^st = nt{t) Xl(-l)''(4n+lC2n+2 + 4n+2C2n+l). (13) 



By solving Eq. ( 11 ) to leading order in Sr, Sb, Se, we find 
that r is given by Eq. (pi) with 

2672 



L Srl32 + 2Sea2 2SeP2 - Sra2 . 



.(14) 
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Here 0^2 = 

A2 = ^{Sb' - Sr' - USef + (2SbSr - USeSb) and 72 
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FIG. 1. Time evolution of r{t) of model (10), in the absence 
of spin orbit (black curve), and for de/dr = 0.3 and F/tq = 0.1 
(red curve). Squares indicate equally spaced time intervals, 
and illustrate the time duration of the transition to r = — 1. 
The loops cross points ±1 at the TRIM 0,T/2, respectively, 
corresponding to a nontrivial Z2 index v — 1. Inset shows the 
time dependence of the ^/-component of e0(t) (black curve) 
and of 0(t) for de/dr = 0.3 and F/tq = 1,0.1,0.01 (blue red 
and green curves, respectively). At F/tq -^ 0, the transition 
becomes sharper and r{t) describes a rotation around a fixed 
axis, e0(T/2), implying that the pump is optimal. 



FIG. 2. Time evolution of the r{t) of model ([l3|, with 
{5e/5T = 0.1, red curve) and without (black curve) spin orbit 
coupling, for F/tq = 0.1. The loops cross two resonances dur- 
ing a cycle and thus belong to the trivial class of pumps. At 
finite spin orbit coupling, r{t) follows two different geodesies, 
with non-parallel tangents at the TRIM, indicated by the 
red arrows. The inset shows the time dependence of the y- 
component of e^(t) (black curve) and of (j){t) for de/dr = 0.1 
and F/to = 0.1. The different values of ecf^iti) at the reso- 
nances, = TT, are marked by the blue circles. 



Fig. [2] shows the time evolution of the reflection ma- 
trix, f(t), for model (13), with (red curve) and without 
(black curve) spin orbit coupling. The loops cross two 
resonances at the TRIM, corresponding to the topologi- 
cally trivial class of pumps, i/ = 0, Eq. ([2|. In contrast to 
the nontrivial class of pumps, these loops are fixed by a 
single point. As a result, they can be contracted to this 
single point by a continuous deformation of the start- 
ing Hamiltonian without violating condition (IT]). This is 
manifested at finite spin orbit coupling, where the loop 
avoids the north pole, f = 1, while remaining fixed to the 
south pole, f = — 1, as illustrated in the red curve. The 
inset shows the time dependence of the ^/-component of 
the rotation axis Ccfyit) and the angle (j){t) for finite spin 
orbit coupling. The angle crosses two resonances, (j) = tt^ 
during a pumping cycle, which correspond to different 
values of the vector e0(ti) (see blue circles in the inset). 
The resulting reflection matrix traces two loops, which in 
the limit F/tq -^ converge to two different geodesies, 
with non-parallel tangents at the TRIM, indicated by 
the arrows, see red curve in main figure. Hence, in the 
presence of a finite spin orbit coupling, the spin pumped 
during a cycle is not quantized, and the pump operates 
with the generation of finite noise. 

Conclusions: — We have studied the class of one- 
dimensional pumps with a time reversal restriction on 
the pumping cycle, ([l]). These systems can be viewed 
as a mapping of two-dimensional time-reversal invariant 
bulk insulators placed on a cylinder, where t corresponds 
to a magnetic flux threading the cylinder. We found that 
the scattering matrices of the pumps are dichotomized by 
the Z2 topological index. We have shown that the class 



of pumps characterized by a nontrivial topological index 
allows, in contrast with its topologically trivial counter- 
part, for the noiseless pumping of quantized spin, even if 
spin is not conserved. This observation sheds light on the 
topological classification of two-dimensional time reversal 
invariant insulators. 
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